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O 

We reproduce the DIS measurements of the proton structure function at high energy from the 
, dipole model in momentum space. To model the dipole-proton forward scattering amplitude, we use 

D ■ the knowledge of asymptotic solutions of the Balitsky-Kovchegov equation, describing high-energy 

QCD in the presence of saturation effects. We compare our results with the previous analysis in 
coordinate space and discuss possible extensions of our approach. 



i. introduction 

On 
O 

■ One of the most intriguing problems in Quantum Chromodynamics (QCD) is the growth of the cross sections for 
t-H , hadronic interactions with energy. As well-known, the increase of energy causes a fast growth of the gluon density 
and consequently of the cross sections. At very high energies, this growth should not continue indefinitely and at 
some point, one has to deal with gluon recombination and multiple scattering in order to restore unitarity. This 
interaction between overlapping partons is called saturation and has deserved active studies over the last thirty years 

JL . More generally, the large amount of work devoted to the description and understanding of perturbative QCD 
in the high-energy limit covers the description of saturation on the theoretical side as well as its applications to 
phenomenology. The theoretical contribution comes mainly from the development of non-linear QCD equations 
describing the evolution of scattering amplitudes towards this limit, together with the search of the solutions to 
those equations. The simplest of such equations is the Balitsky-Kovchegov (BK) equation 0, Q, which corresponds 
to the Balitsky-Fadin-Kuraev-Lipatov (BFKL) Q linear evolution equation with the addition of a non-linear term 
responsible for the saturation of the growth of gluon density. It has been shown [10( that the BK equation is in the 
equivalence class of the Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) nonlinear partial differential equation [ll[ , 
which admits travelling-wave solutions, translating, in terms of QCD variables, into geometric scaling as we shall 
explain below. 

From the phenomcnological side, the geometric scaling has been observed at the DESY ep collider HERA, in the 
measurements on inclusive j*p scattering 12]. This phenomenological feature of high-energy Deep Inelastic Scattering 
(DIS) is expressed as a scaling property of the virtual photon-proton cross section 



- r— I 

X 



that is, the cross section depends on the scaling variable r = Q 2 /Q 2 (Y) instead of Q 2 and Q 2 (Y) separately. Here 
Q 2 is the virtuality of the photon, Y = log(l/x) is the total rapidity, x is the Bjorken-x, related to the centre-of-mass 
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energy through s = Q 2 x and Q S (Y) is an increasing function of Y called the saturation scale. The geometric scaling 
is actually equivalent to the formation of travelling-wave solutions for the BK equation. This is thus a remarkable 
consequence of saturation, which extends arbitrarily far beyond the fully saturated domain, i.e. in the dilute regime 
where saturation effects may seem negligible. 

In this paper, we use the dipole model [13j to relate the 7*p cross-section to the dipole-proton forward scattering. 
This approach has already been proven successful e.g. in [bH. flBL [HI] . Our approach here is to parametrise the dipole- 
proton amplitude in momentum space, where the travelling-waves have been originally derived. We shall discuss the 
advantages of our method and compare it with previous results in the literature later on. 

The plan of this paper is as follows. In section [Til we relate the 7*p cross section to the dipole-proton scattering 
amplitude within the dipole framework. We then discuss, in Section IIIII how one can describe the dipole scattering 
amplitude from the properties of the BK equation. In Section I IV Al we gather all information to build the complete 
model for the proton structure function. The fitting procedure used to compare our model with the experimental 
measurements is explained in Section llVBI and the results of the fit are presented in Section ftV CI We discuss the link 
with previous approaches in the literature and possible situations in which our work can find interesting applications 
in Section Ivl 

II. THE DIPOLE MODEL 




FIG. 1: Picture representing the dipole model. 



We consider the collision between a virtual photon and a proton at high energy. In a frame where the photon 
travels fast, one can consider that it fluctuates into a qq dipole. The lifetime of this dipole being much longer than 
the time of interaction with the proton, one can write the cross section as a product of the wavefunction for a photon 
to go into a dipole times the dipole-proton cross section, as shown on figure [TJ It leads to the well-known formula 

4l(Q 2 X) = / d 2 r j\z \^ T , L (r,z;Q 2 )\ 2 a^(r,Y), (2) 

where crj-^ir, Y) is the dipole-proton cross section. The transverse and longitudinal photon wavefuctions in this 
expression are computable in perturbative QED. They are given by 

|* T (r, z; Q 2 )\ 2 = e ? { I** + C 1 " *)"] W(Q 9 r) + <K 2 {Q q r)} (3) 

9 

\* L (r,z;Q 2 )\ 2 = ^i^e 2 {4QV(l-z) 2 if 2 (g g r)} (4) 

q 

where N c is the number of colours, Q 2 = z(l — z)Q 2 + m 2 , m q the mass of the quark of flavour q and iio,i are the 
Mc Donald functions of rank zero and one, respectively. This expression extends the Kovchegov formula (Eqs. (4), 
(5a) and (5b) in the second reference of [8|) by including the mass of the quarks. 

If one treats the proton as a homogeneous disk of radius R p , the dipole-proton cross section in eq.© is usually 
taken to be proportional to the dipole-proton forward scattering amplitude T(r, Y) through the relation 



o-£ p p (r,Y) = 2nR 2 p T(r,Y). 
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The proton structure function F% can be obtained from the j*p cross section through the relation 

12 

- ''•/■ 



F 2 (x,Q 2 ) 



4tt 2 



4*>,Q 2 ) + af p (:r,Q 2 ) 



(5) 



As we shall see in the next section, the BK equation describes the high-energy evolution of the dipole-proton scattering 
amplitude T. However, the asymptotic behaviour of its solutions is naturally expressed in momentum space. Hence, 
we want to express the 7*p cross section in terms of T(£:, Y), the Fourier transform of T(r, Y): 



1 f d 2 r f°° dr 

T(k, Y) = — / — e 4kr T(r, Y) = / -J (kr) T(r, Y). 
2n J r z Jo r 

After a bit of algebra, we obtain the following expression relating the proton structure function to T(k,Y): 



F 2 {x 1 Q 2 ) 



2 _ Q 2 R 2 p N c dk f 1 



4tt 2 







dz|*(fc 2 ,z;Q 2 )| 2 T(fc,F), 







where the wavefunction is now expressed in momentum space 
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Having described the DIS structure function through l[7|). the approach used to describe the scattering amplitude 
T(k, Y), the forward scattering amplitude in momentum space, will be presented and discussed in the next section. 



III. SCATTERING AMPLITUDES IN HIGH-ENERGY QCD 

Let us now consider a fast-moving colourless qq dipole of transverse size r = |x— y|, where x and y are the coordinates 
of the quark and antiquark respectively, interacting with a given dense target. In the large- N c approximation (N c 
being the number of colours), and in the mean-field approximation, the high-energy behaviour of the dipole forward 
scattering amplitude T(x, y;Y) follows the BK equation 0, In coordinate space this equation reads 

dyT(x, y; Y) = [ d 2 z { * ~ y)2 [T(x, z; Y) + T(z, y; Y) - T(x, y; Y) - T(x, z; Y)T(z, y; Y)] (9) 
2tt J (x - z)-(z - y) 2 

where a — a s N c /n, a s is the strong coupling constant, considered fixed. If one neglects the dependence on the impact 
parameter b = (x + y)/2 and integrates out the remaining angular dependence of r, ([9]) becomes an equation for 
T{r, Y). The latter can be expressed in momentum space using ([6]). 
One finds T(k, Y) obeys the BK equation in momentum space 

d Y T = a X {-d L )T ~aT 2 , (10) 

where a = a s N c /ir and 

X(7) = 2^(1) - V>(7) - 1>(1 ~ 7) (11) 

is the characteristic function of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) kernel 0] and we used L — log(fc 2 /fc§) 
with k some fixed soft scale. 

In our approach, we do not want to use directly numerical solutions of (|10p as the input for the dipole scattering 
amplitude but rather use our knowledge of the properties of its solutions to build an analytical expression for T{k, Y). 
As we shall explain in more details in the next section, this allows, for example, to take into account next-to-leading 
order corrections to the BFKL kernel, which provides a much better description of the data. 

It has been shown [l(| that, after the change of variables 



t = aY, x~\og(k 2 /kl), u(x,t) oc T(k,Y), 



(12) 
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the BK equation reduces to the Fisher-Kolmogorov-Petrovsky-Piscounov (F-KPP) equation [1JJ for u when its kernel 
pip is approximated in the saddle point approximation i.e. to second order in the derivative d^, the so-called diffusive 
approximation. The F-KPP equation is a well-known equation in non-equilibrium statistical physics, whose dynamics 
is used to describe many reaction-diffusion systems in the mean-field approximation: 



d t u(x, t) — d 2 u(x, t) + u(x, t) — u 2 (x, t), 



(13) 



where t and x are, respectively, the time and space variables. 

This equation has been extensively studied since more than sixty years and, in particular, it is known to admit 
travelling waves as asymptotic solutions. This means that the solution u(x, t) to equation Q13[) takes the form u(x~v c t) 
of a front travelling to large values of x at a speed v c without deformation. 

This property of the F-KPP equation is actually true if one considers the BK equation (flU)) with the full BFKL 
kernel [1 01 ] - At asymptotic rapidities, the amplitude T(k, Y), instead of depending separately on k and Y , depends only 
on the scaling variable k 2 /Q 2 {Y), where we have introduced the saturation scale Q 2 (Y) = k 2 exp(v c Y), measuring 
the position of the wavefront (more precisely, it is L S (Y) = log^ 2 /^ 2 .) = v cY which measures its position). 

A more detailed calculation allows also for the extraction of two additional subleading corrections, resulting into 
the following expression for the tail of the scattering amplitude 



T(k,Y) 



Q 2 s(Y) 



log 



Q 2 s(Y) 



exp 



log 2 (k 2 /Q 2 s (Y)) 
2a X "(lc)Y 



(14) 



where x" denotes the second derivative of the BFKL kernel (jlip with respect to 7 and the saturation scale 



Q 2 S (Y) =fc 2 exp av c Y- 



^ logon -A 

27c 7c 
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(15) 



The critical parameters j c and v c are obtained from the knowledge of the BFKL kernel alone and correspond to 
the selection of the slowest possible wave: 



v c — mma- 

7 



.x(t) -X(7c) 



7 



7c 



«x'(7c)- 



= 4. 



(16) 



For the leading-order BFKL kernel (TTTj) . one finds 7 C = 0.6275..., and v c 

Physically, the property of geometric scaling has deep consequences. It expresses the fact that when one moves 
along the saturation line, the behaviour of the scattering amplitudes remains unchanged. In addition, saturation 
introduces naturally the scale Q s which provides an infrared cut-off solving the infrared instability problem of the 
BFKL equation. The subleading corrections in (T4"|) also have a very important role. The last term introduces an 
explicit dependence on the rapidity Y and hence violates geometric scaling. However, this term can be neglected 
when 



log 2 (k 2 /Q 2 (Y)) 



2a X "hc)Y 



< 1. 



This means that geometric scaling is obtained for 



log(fc 2 /Q 2 (F)) < ^2 X "( 7c )ar, 



i.e., in a window which extends \JY above the saturation scale. It is a remarkable property that, at high-energy, the 
consequences of saturation are observed arbitrarily far in the tail, where T(k,Y) is much smaller than 1. 



IV. DESCRIPTION OF THE 7 *p DATA 

The method that will be used to describe the experimental measurements of the proton structure function can be 
split in three steps. We shall first build a QCD-based model for the scattering amplitude, which can directly be used 
in ([7| to obtain F2. Then the details concerning the remaining ingredients needed to obtain F2 are given and we 
specify the dataset. We finally present the results. 
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A. Dipole scattering amplitude 



Expression (|T4|) only gives a description of the tail of the wavefront T(k,Y) -C 1 (k 3> Q s ). In order to complete 
the description, we also need expressions for T around the saturation scale and at saturation. In the infrared domain, 
one can show, for example by computing the Fourier transform @ of a Heaviside function T{r) = Q(rQ s — 1), that 
the amplitude behaves like 



T 



Q.{YY 



Y\ k< ^c- 



log 



Qs(Y) 



(17) 



where c is an unfixed constant. 

Those parametrisations (Equations (fT4|) and (fT7|) ) describe fully the asymptotic behaviour of the amplitude and we 
are left with the matching around the saturation scale. The easiest way to do that is probably to use (|14| for k > Q s 
and (TIT)) for k < Q s and match the constant c to obtain a continuous distribution. However, this definition by parts 
would certainly introduce oscillations in the coordinate space T(r, Y) which may even lead to negative amplitudes. 
Thus, the best way to obtain the description of the transition to the saturation region is to perform an analytic 
interpolation between both asymptotic behaviours. 

In order to obtain an interpolation model which describes the transition from the dilute regime to the saturation 
one, the idea is to build the latter domain from the former. Our starting point is an expression which is monotonically 
decreasing with L and which reproduces (up to the logarithmic factor), the amplitude for geometric scaling (114 



T di i = exp 



-7c log 



A- 2 



Q 2 s(Y) 



4 2 ed-l0g 2 (2) 



2a X "{lc)Y 



(18) 



with 



i rod = log 



and Qi(Y) = k£e 



2 &v a Y 



(19) 



This result is unitarised a la Glauber-Mueller i.e. Turrit = 1 — exp(— Tdii) and we finally need to reinsert both 
logarithmic behaviours in the infrared and in the ultraviolet. We obtain that the following choice gives good results: 



T(k,Y) 



log 



Qs 
k 



l 



l-e" 



(20) 



Equations {T5}, l|19p and (|2"D|) determine our model for the scattering amplitude. They provide a valid and simple 
interpolation between the QCD constraints (|14p and l|17p. They will be considered in order to describe the structure 
function F2 given by 



B. Dataset 

We fit all the last HERA measurements of the proton structure function from HI [TU, ZEUS [lj|, with our analysis 
being restricted to the following kinematic range: 

[x < 0.01, 

1 0.045 <Q 2 < 150 GeV 2 . v ; 

The first limit comes from the fact that our approach is meant to describe the high-energy amplitudes i.e. the small x 
behaviour. The second cut prevents to reach too high values of Q 2 for which DGLAP corrections need to be included 
properly. This gives a total amount of 279 data points. In addition, we have allowed for a 5% renormalisation 
uncertainty on the HI data. 

Concerning the parameters, we have kept 7 C = 0.6275 and a = 0.2 fixed although a particular choice for a s only 
results in a renormalisation of the other parameters. For the convolution with the photon wavefunction, we have 
assumed different situations for the quarks masses: the light-quarks mass m q has been set to 50 or 140 MeV while, 
we have used m c — m q or m c = 1.3 GeV for the charm mass. This leaves v c , x", k 2 , and R p as free parameters. 
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FIG. 2: Results for the fit to the virtual photon-proton cross section. The data points and the fit are shown as function of Q 
for different values of x. Those x values are indicated for each curve. For both plots, the successive curves have been rescaled 
by powers of 2 (1 to 2~ 9 from top to bottom) to ensure clarity. The plus corresponds to the ZEUS points, the cross to HI. 



Masses 


kl (10" 3 GeV 2 ) 


v c 


// 

Xc 


R p (GeV) 


X 2 /nop 


m, = 50 MeV, m c = 50 MeV 
m, = 50 MeV, m c = 1.3 GeV 
m, = 140 MeV, m c = 1.3 GeV 


3.782 ±0.293 
7.155 ±0.624 
3.917 ±0.577 


1.065 ±0.018 
0.965 ±0.017 
0.807 ± 0.025 


4.691 ±0.221 
2.196 ±0.161 
2.960 ± 0.279 


2.770 ± 0.045 
3.215 ±0.065 
4.142 ± 0.167 


0.960 
0.988 
1.071 



TABLE I: Results from the fit to the F2 data. The values of the parameters with their respective errors are indicated, together 
with the x 2 P er data point. 



C. Results 



The parameters obtained from the fit are shown together with the x 2 per point in Table Q] We have also plotted 
the comparison with the experimental data for F% on figure [2] as well as the charm structure function i 7 ^ on figure 
[3l For both cases, the curve correspond to m q — 50 MeV and m c = 1.3 GeV. We verify a good agreement with the 
measurements of F% due to the small \ 2 provided by the fit. Moreover, the F% predicted by the parametrisation is in 
reasonable agreement with the experimental results, which shows the robustness of the model proposed in this work. 

Within our parametrisation, the saturation scale Q s corresponds to the energy-dependent scale at which the dipole 
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FIG. 3: Predictions for the Hl[2(j and ZEUS[2l|] measurements of the charm structure function. The result is presented as a 
function of x for different values of Q 2 . For each subplot, the value of Q 2 is given in GeV 2 . 

scattering amplitude 1 is T(k = Q S {Y), Y) = [l + log(2)](l - l/e) » 1.07. For the fit corresponding to m q = 50 MeV 
and m c — 1.3 GeV, we obtain a saturation scale Q s = 0.206 GeV for x = 10~ 4 and Q s = 0.257 GeV for x = 10 -5 . 
Although these values may seem rather small, we have to emphasise that they correspond to large values for T. If, 
instead, we extract the saturation scale by requiring that T — 1/2 when k = Q s , we get Q s = 0.296 (resp. Q s — 0.375) 
GeV for x = 10~ 4 (resp. x = 10~ 5 ). These last values are still a bit smaller than the saturation scales observed in 
previous studies (Q s f» 1 GeV for x w 10 -5 ) and, hence, tend to confirm the tendency for the saturation scale to 
decrease when heavy-quark effects are taken into account. 



V. CONCLUSIONS AND DISCUSSIONS 

In this work we have investigated the travelling-wave solutions of the BK equation which describe the forward 
scattering amplitude at high energies and tested their phenomenological implications for the virtual photon-proton 
scattering. We have proposed an expression for the amplitude in momentum space (see eqs. (|18l) . (|19[) and (|20[1 ) which 
interpolates between the behaviour of the dipole-proton amplitude at saturation and the travelling- wave, ultraviolet, 
amplitudes predicted by perturbative QCD from the BK equation. This expression was used to compute the proton 



Note that, when working in momentum space, the logarithmic behaviour of the amplitude in the infrared allows it to take values larger 
than 1. This is not a contradiction with unitarity. 
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Parametrisation 


Quark masses 


nop 


X 2 /nop 


GBW [14] 


m, = 140 MeV, m c = 140 MeV 


372 


1.5 


IIM [15] 


m q — 140 MeV, no charm 


156 


0.81 


BGK [16] 


m q = MeV, m c = 1.3 GeV 


288 


1.06 


This work 


m q = 50 MeV, m c = 1.3 GeV 


279 


0.988 



TABLE II: Comparison between the results from the parametrisations in [bH . HH Il6j and the the present work. 

structure function F2 (in the framework of the dipole model) and tested against the HERA data. We obtained a good 
fit with light quark masses m u ^d, s = 50, 140 MeV and heavy charm mass m c = 1.3 GeV. 

At this point, it is interesting to compare our results with those from previous approaches in the literature. In the 
last decade some dipole models of DIS at small- a; have been proposed and turned out to be successful in describing 
the experimental data. Among these, our comparisons will be particularly focused on the pioneering Golec-Biernat- 
Wusthoff (GBW) model [ijj], the lancu-Itakura-Munier (IIM) model [l5| and recent developments concerning the 
Bartels-Golec-Biernat-Kowalski (BGK) model [HI, which consists of the GBW model improved by the incorpora- 
tion of a proper gluon density evolving according to the DGLAP evolution equation. Some results from these models 
are shown in Table HT1 For each parametrisation, we indicate the values used for the (light and charm) quark masses, 
the number of points and the respective values of x 2 P er data point obtained from the fit. 

First of all, it should be stressed that these models were developed in coordinate space and not in momentum space, 
as in the present work. Our choice is directly motivated by the analysis of the BK equation in momentum space leading 
to universal asymptotic results on which we heavily rely. Although most approaches are able to reproduce the F% 
data with a good x 2 , our model can be differentiated from the previous ones at two levels. On the one hand, our 
analysis is based on the BK equation to account for unitarity effects. Thus, we expect it to be more precise than the 
Glauber-Mueller- like approach used in [3] and [HI , especially, in the small- a; and low Q 2 domain under study. On the 
other hand, we improved the IIM model by including massive charm. Moreover, in both GBW and IIM models, one 
finds some difficulties after performing their Fourier transform, that is, changing from the coordinate to momentum 
space. Such aspects were analysed in details in (22J. In the case of GBW model, one obtains a Fourier transform of 
the dipole cross section which presents an unrealistic perturbative behaviour, while in the case of IIM it presents non- 
positivity values. These problems are tamed in our model, where the inverse Fourier transform (scattering amplitude 
in coordinate space) remains between and 1. The resulting dipole cross section presents the colour transparency 
property, i.e. a dip ~ r 2 when r — > and the saturation property, i.e. a dip ~ at large r. 

These features associated with the good data description and small x 2 provide that the dipole scattering amplitude 
proposed in this work should be a good parametrisation to investigate the properties of the observables at RHIC and 
LHC energies, considering the dipole approach. 
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